Creation of two vortex-entangled beams in a vortex beam collision with a plane wave 
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Physics of photons and electrons carrying orbital angular momentum (OAM) is an exciting field 
of research in quantum optics and electron microscopy. Usually, one considers propagation of these 
vortex beams in a medium or external fields and their absorption or scattering on fixed targets. Here 
we consider instead a beam-beam collision. We show that elastic scattering of a Bessel vortex beam 
with a counterpropagating plane wave naturally leads to two vortex-entangled outgoing beams. 
The vortex entanglement implies that the two final particles are entangled not only in their orbital 
helicities but also in opening angles of their momentum cones. Our results are driven by kinematics of 
vortex-beam scattering and apply to particle pairs of any nature: ey, e + e~ , ep, etc. This collisional 
vortex entanglement can be used to create pairs of OAM-entangled particles of different nature, and 
to transfer a phase vortex, for example, from low-energy electrons to high-energy protons. 



PACS numbers: 42.50.Tx, 41.85.-p, 03.65.Nk, 12.20.-m 

Introduction. — Laser beams carrying non-zero orbital 
angular momentum (OAM) are well-known and routinely 
used in optics, [TJ [2] - Wavefronts of such a beam are not 
planes but helices, and each photon in this light field (a 
twisted photon) carries a well-defined OAM quantized in 
units of h. Applications of twisted photons range from 
microscopy to astrophysics, [5] . They are also interesting 
for quantum information science because of the high di- 
mensionality of the OAM state space and the possibility 
to create OAM-entangled pairs of photons, [HH]. 

Wave-fronts with phase vortices can exist for electrons 
and other particles as well. Recently, following the sug- 
gestion of [6], electron beams carrying OAM were exper- 
imentally demonstrated, first using phase plates [7j and 
then with fork diffraction holograms [8]. Such electrons 
carried kinetic energy as high as 300 keV and the orbital 
quantum number up to m ~ 100. Exact solution of the 
Dirac equation representing a relativistic electron vortex 
was recently derived in [5]. 

Experimental realization of vortex beams opens up 
the possibility to study head-on beam-beam scattering, 
with one or both colliding beams carrying OAM. Us- 
ing sufficiently energetic twisted electrons and photons 
in &y, ee, e + e _ , ep or electron-nucleus collisions, one can 
probe quantum electrodynamics and perhaps hadronic 
processes in a novel way, [10]. In addition, the Comp- 
ton backscattering of optical twisted photons from an 
ultrarelativistic electron beam was suggested in [TT] to 
generate twisted photons in the GeV energy range. 

All such processes share certain universal features, 
which are driven by kinematics of free-space beam-beam 
collision and are not sensitive to its microscopic nature. 
Some of them were studied in [TO] [12]. One particular 
finding was the necessity of the use of orbital helicity, 
the OAM projection on the average propagation direc- 
tion (which is different for each particle) rather than on 
a fixed common axis. This notion becomes especially im- 
portant for high-energy and for high-angle (strongly non- 
collinear) scattering, as it removes strong "instrumental 



OAM" effects caused by an unfortunate choice of the ref- 
erence axis and leaves us with the "intrinsic OAM" of 
the twisted state. The usage of orbital helicity (without 
calling it so) was also advocated in [T3]. In the present 
paper whenever we refer to OAM we actually mean the 
orbital helicity. 

In this Letter we explore another universal feature of 
vortex beam collisions. We study whether the outgoing 
wavefronts contain vortices and what are the correspond- 
ing orbital helicities mi and m^- This issue has never 
been addressed in a generic set-up; in previous calcula- 
tions of vortex beam scattering at least one of the final 
particles was assumed to be a plane wave. We find that 
the final particles are vortex- entangled: that is, they are 
entangled both in the (mi, 777,2 )-space and in the space of 
opening angles. Since for non-collinear scattering there 
is no total orbital helicity conservation law, selecting one 
specific mi does not lead to a unique mi- However by 
using the other degree of freedom, the opening angle, one 
can collapse the second particle to a more-or-less definite 
m2, with important practical implications. 

Note that similar issues arise in the context of OAM- 
entangled photon pairs production by spontaneous para- 
metric down-conversion (SPDC) in a non-linear medium, 
see for example [131 [U]. In the particle language, such 
medium induces decay of the pump photon into two pho- 
tons of lower energy. Although a convenient source of 
OAM-entangled photon pairs, this process is specific to 
photons and its dynamics is defined (and limited) by the 
exact properties of the crystal slab, including its thick- 
ness. In contrast to that, the entanglement discussed in 
this Letter is universal (not specific to particle species), 
it spontaneously occurs in a free-space collision (and not 
induced by a medium), and it allows one to study its 
behavior in the entire available kinematical range. 

Describing twisted states. — To simplify the calcula- 
tions, we describe the twisted states as paraxial Bessel 
beams. We consider the scalar case only (polarization 
degrees of freedom can be straightforwardly incorporated 
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in the paraxial approximation) and use the conventions 
of [11]. All transverse vectors are given in bold, while 
three-vectors will be indicated by the arrow. A Bessel 
twisted state is a non-plane wave solution of the free 
wave equation with a definite frequency uj, longitudi- 
nal momentum k z , modulus of the transverse momen- 
tum |k| = x and a definite z-projection of orbital an- 
gular momentum m. When written in cylindric coor- 
dinates r,ip r ,z, it has form |x, m) = e~' lu}t+ ' lk ' z ip ><:m (r), 
with VVm( r ) = e" nVr J m (xr ) y/x/2-K, where J m (x) is the 
Bessel function. A twisted state can be represented as a 
superposition of plane waves 



\x, m) = e 



d 2 k 



cw(k)e lkr , (1) 



where 



<w(k) = (-i) m e m ^V^ ' 5(|k| / - X) . (2) 

Thus, the allowed momenta lie on the edge of a cone 
with the opening angle arctan(x/£; z ). More properties of 
twisted states, including their normalization and phase 
space density, can be found in jTU], [TT]. Here we only 
note that Eq. ([I]) in fact describes the passage from plane 
waves to twisted particles in description of a scattering 
process. 

A Bessel state |x, m) with fixed x is non-normalizable 
in the transverse plane. A much better approximation 
to physically realizable states such as Bessel-Gaussian or 
aperture-limited beams is given by a fixed-m superposi- 
tion of Bessel states 



\x ,a;m) = / dx/(x)|x,m) 



(3) 



with a properly normalized weight function /(x) peaked 
at xq and having width a. This state is normalizable 
(and localized) in the transverse plane and is assumed 
to be monochromatic (k z is supposed to vary with x so 
that the energy is constant). Properties of such states 
and their important role in resolving the non-forward-to- 
forward paradox in Bessel beam scattering was discussed 
in PI] . 

Vortex beam scattering. — To describe scattering of the 
Bessel state \k, m) and a plane- wave, we start with a 
generic elastic two-particle scattering in the usual plane 
wave basis. The initial particle four-momenta are de- 
noted as k and p, the final four-momenta are k\ and k%. 
The scattering matrix element of this process is repre- 
sented as 



S 



PW 



i{2ir) 4 6( 4) {k+p-k 1 -k 2 )M, 



(4) 



where the amplitude A4 is calculated according to the 
standard Feynman rules. Then we pass from the plane- 
wave to twisted scattering by convoluting Q with 
the weight function o Ximj , for each initial Bessel- 
beam twisted state, and with a%. m . for each final 



twisted state. For example, if only one initial parti- 
cle is twisted, the scattering matrix element is St w = 
J d 2 ka >C m(k)Spw /(2n) 2 . Owing to the delta-function 
inside a^ m , this representation contains only one integra- 
tion with respect to the azimuthal angle tp of the trans- 
verse momentum k. This integration is eliminated by the 
transverse momentum delta- function in Spw, Q- If we 
assume for simplicity that p = (0,0, p z ) and introduce 
the well-defined transverse vector k 12 = kx + k 2 with 
modulus ki2 and azimuthal angle <fi 2 , then 



Stw oc 



d 2 k 



(2<p 

(-0* , 

(2tt) 3 / 2 



a^ m (k)j( 2 )(k-k 1 -k 2 )M(k) 



(5) 



Note that the scattering amplitude M. itself is not inte- 
grated but is just taken at a specific value of the initial 
momentum. 

Do the wavefronts of the two outgoing particles con- 
tain phase vortices? Eq. ^ provides the answer to 
this question. Let us first suppose that the final mo- 
mentum k 2 with modulus |k 2 | and azimuthal angle ip% 
is fixed (the second particle is projected on a plane 
wave). Then, the final momentum of the first parti- 
cle ki is not uniquely defined but belongs to a circle of 
radius x around the point — k 2 . Its modulus changes 
in the interval ||k 2 | — x\ < |ki| < |k 2 | + x, and for 
a given |ki| its azimuthal angle tpi takes two values: 

k 2 -k 2 \ 

The azimuthal angle 
also 



if2 ± arccos 



/ ~ 2 -k 2 -k 2 \ 

V 2|kx||k 2 | )■ 



if I — y-z -l. ^ 2|ki||k 2 | 

ipi2, which parametrizes the points on the circle 



(2£ 2 +k 2 — k 2 \ 
- 2>c |k 2 | 1 ) ' anc ^ ^ ne 

plus-minus signs in these two expressions are correlated. 
When |ki| spans the allowed interval, the angle <pi2 cov- 
ers the entire circle, although the intensity of the scat- 
tering is modulated by the amplitude Ai. 

Note that the scattering amplitude for k x = — k 2 is 
exactly zero for any value of x. Passing to the three- 
dimensional vectors, one can define the axis fti\\ (k) +p — 
fc 2 , where (k) = (x,m\k\x,m), and claim that scatter- 
ing exactly in this direction in absent even of the initial 
twisted state is represented by a transversely localized 
state ((3j). Therefore, n\ represents the direction of the 
phase vortex (the line of zero intensity and undefined 
phase). It is, therefore, natural to expect that at fixed fc 2 
the final state of the first particle can be approximated 
by the twisted state |xi, mi) defined with respect to this 
direction, with xi ~ x and with orbital helicity m\ ~ m. 
For the specific kinematics of the Compton backscatter- 
ing, this expectation was confirmed in [12] . 

It is clear that if one fixed the final momentum k± of the 
first particle, one would rederive similar conclusions for 
fc 2 : the outgoing wave of the second final particle would 
contain a phase vortex in the direction of n 2 ||(/c) + p — 
k±. Therefore, there can be no unambiguous assignment 
which of the final particles is twisted and which is not. 
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FIG. 1: (Color online.) Kinematics of the triple-twisted scat- 
tering in the c.m.s. 

Triple-twisted scattering. — The "plane wave/twisted 
state" entanglement derived above rests on the assump- 
tion that one of the final particles is a plane wave. Below 
we show that, in general, the final state in this elastic 
collision can be represented by an entangled state of two 
vortex beams. The final twisted particles are entangled 
not only through m's, but also through x's, and the pre- 
viously found "plane wave/twisted state" entanglement 
is just a particular case of this general result. 

Let us again consider the elastic collision of a Bessel 
state | x, m) with the plane wave with momentum p and 
represent the two final particles as Bessel states \n\,m\) 
and \x 2 ,m 2 ). Kinematical conventions are shown in 
Fig. [T] In order to simplify the calculations, we con- 
sider this process in the c.m.s. frame, where (k) + p = 0, 



Here we introduced parameter £ via sin£ = q/x, |£| < 9, 
and the following angles 




Besides, A is the area of the triangle with sides x\, 
x 2 , and x = >rcos£, and Si and S 2 are two of its in- 
ner angles: cosJi.2 = (x 2 + x\ 2 — x\ 1 )/2xxi i 2- The 
results of [T^] for the case when one of the final parti- 
cles is a plane wave can be recovered from ^ in the 
case of m 2 = and x 2 — > 0, using the following re- 
lations, [in]: |PW(k 2 )) = lim^o y/2ir/x 2 \x 2 ,0) and 
lim K2 _y A -1 = 4ttS(x 2 - h\). 

Our main result Q allows us to characterize the en- 
tanglement of the two final particles. They are entangled 
not only in the (mi,m 2 ) space, but also in the (xi,x 2 ) 
space. Existence of the triangle with sides xi, x 2 , and x, 
means that the allowed values of (xi,x 2 ) lie inside the 



and choose the common quantization axis z', character- 
ized by the polar angle 9, for both final particles. The 
two axes, z and z', define the scattering plane; working 
in this plane, we choose axis x as the one orthogonal 
to z, and axis x' as orthogonal to z' . Without loss of 
generality, we assume that the azimuthal angle ip for the 
initial twisted state is measured from axis x, while the 
azimuthal angles ip± and ip 2 of the final twisted states 
(which lie in the plane orthogonal to z') are measured 
from the axis x' . The average values of the final mo- 
menta (ki) and (k 2 ) are parallel to z'\ we denote their 
difference as q = |(fci)| — \ (k 2 }\ = k\ z i + k 2z >. 

The scattering matrix element for such a triple-twisted 
scattering, S3tw, is equal to 

f d 2 k d 2 k x d 2 k 2 

J (27T)2 (27T)2 (2tt)2 MKj^m^lja^ (k 2 )b PW . 

(6) 

The three-dimensional delta-function 5^ (k + p — k\ — k 2 ) 
inside Spw eliminates all three integrals over the az- 
imuthal angles (f, ifi and (p 2 setting them to certain val- 
ues. The amplitude M. does not affect the integral and 
is just taken at the momenta fc, k\ and k 2 correspond- 
ing to these azimuthal angles. In order to show the key 
features of the results, we now assume that the ampli- 
tude M. is a slowly varying function, and in the paraxial 
approximation we take Ai(ki) s» A4((ki)) = Aio out of 
the integral. Then, the integral can be computed exactly, 
and the triple-twisted scattering matrix element equals 



(7) 

I 

stripe defined by 

\x\ — x 2 \ < x < x\ + x 2 . (9) 

As for the m\- and 7712-distributions, they are infinitely 
wide for pure Bessel states. However if the initial and 
final twisted states are taken as transversely localized 
states Q with peak positions x oi and widths cr^, then the 
orbital helicity distributions become rather narrow with a 
typical width of 0{xQi/ ai). Their shapes also depend on 
the location of the central point (xoi, x$ 2 ) in the stripe 
of allowed values ([9| . For illustration, we show in Fig. [2] 
a typical (m 1; m2)-distribution of the scattering intensity 
as a proportional box plot at 9 = 0.2 w 11° for to = 5 
and for an asymmetric choice of >r's: x$ = xqi = 2xq 2 
(the overall scale is arbitrary), er, = xoi/5. This scat- 
tering intensity is obtained by squaring the /^-weighted 
S'-matrix element ([7| and integrating the result over the 
allowed region of q. The (to 1; TO 2 )-correlations seen in 
Fig. [2] arise due to the orbital helicity entanglement of 



J3tw 



MEf-Ei) 



■m 1 +m2 — m 



2 I x\x 2 cos [my;* — (m,\ — m 2 )ip*] ■ cos[mi5i + m2<52] 
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A 



y/ sin 2 9 - sin 2 \ 
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the two vortex beams. 
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FIG. 2: Relative intensity of scattering on the (mi,m2)-plot 
integrated over q for kq = xoi = 2><o2, Oi = Xi/5, m — 5, and 
scattering angle 6 = 0.2. 

The plot also illustrates the absence of the strict total 
orbital helicity conservation rule, although it possesses 
a sizable "ridge" near mi — mi — m. Another feature 
which is visible in the plot is that the final particle with 
the larger x tends to carry larger orbital helicity. This 
property can be read off Q. If < « x, then 
5i <C 1, while 62 can be large. Therefore, mi can be 
large, while any non-zero 7712 will be suppressed by the 
second cosine after integration over certain Xi domains. 
In the limit of vanishing X2 or x% one recovers the "plane 
wave/twisted state" entanglement discussed above. 

Discussion. — Our results bear several implications. 
First, we showed that elastic scattering of a vortex beam 
with a plane wave spontaneously leads to two vortex- 
entangled beams. This can be used to create OAM- 
entangled pairs of particles of different nature (electro- 
photon, electron-proton, etc.), including those for which 
vortex beams are not yet available. So far, only photon 
pairs entangled through OAM were created, [4j [5] . Ex- 
perimental study of this entanglement is feasible with to- 
day technology, for example, by colliding electrons from 
two electron microscopes in the common focal spot and 
detecting the orbital helicity of scattered electrons. Cur- 
rent discussion about OAM conservation during SPDC 
[T3j might also benefit from these experiments. 

Second, by colliding a vortex electron beam with a 
high-energy proton beam with sufficient transverse coher- 
ence and by filtering only scattered electrons with small 
x, one can create energetic twisted protons. This would 
be a unique possibility to generate high-energy particles 
carrying OAM, because fork diffraction grating or phase 
plates are of no use for high-energy (GeV range) parti- 
cles. Due to strong momentum imbalance, the kinemat- 
ics of this process in the laboratory frame differs from 
the c.m.s. example considered here. However for small 
proton scattering angles the integral |6| can also be cal- 



culated analytically, and the OAM-entanglement pattern 
can be studied. 

Third, if a vortex electron beam experiences multiple 
scattering off stray atoms or other particles during its 
propagation, its twisted state can deteriorate. From the 
calculational point of view this process is different from 
propagation of twisted light through a turbulent medium 
or from motion of twisted electrons through chaotic elec- 
tromagnetic fields. In contrast to scattering in an ex- 
ternal potential, the scattering matrix element for a two 
plane wave collision contains the three-dimensional delta- 
function describing momentum conservation, Q. There- 
fore, projecting the initial and final particles on twisted 
states will lead to different results in these two cases. 

In summary, we found that generic elastic scattering 
of a vortex beam with a plane wave leads to two vortex- 
entangled outgoing beams. This feature was ignored in 
all previous calculations of vortex beam scattering, [101 — 
ITS] . The final particles are naturally entangled both 
through their orbital hclicities rrii and through Xi. These 
results are driven by kinematics of vortex-beam scatter- 
ing and apply to particles of any nature (cy, e + e~, ep, 
etc) . They can be used to create pairs of OAM-entangled 
particles of different nature, to transfer a phase vortex to 
high-energy particles (for example, protons), and are im- 
portant for manipulating the electron vortex beams. 
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